Tunneling into a Luttinger liquid revisited 
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We study how electron-electron interactions renormalize tunneling into a Luttinger liquid beyond 
the lowest order of perturbation in the tunneling amplitude. We find that the conventional fixed 
point has a finite basin of attraction only in the point contact model, but a finite size of the contact 
makes it generically unstable to the tunneling-induced breakup of the liquid into two independent 
parts. In the course of renormalization to the nonperturbative-in-tunneling fixed point, the tunneling 
conductance may show a nonmonotonic behavior with temperature or bias voltage. 
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Electron tunneling into a correlated many-electron sys- 
tem is one of the most essential tools to probe the na- 
ture of the correlations. A key concept here is that the 
tunneling density of states reflects how difficult it is for 
electronic states to rearrange themselves to accommodate 
the extra charge of the tunneling electron. The hallmark 
of strong correlations is the "zero-bias anomaly" (ZBA) 
[l[ — the nonlinear behavior of the tunneling current as a 
function of the bias voltage — resulting from a singularity 
in the tunneling density of states at the Fermi energy. 

The prototype model and the best-understood example 
of a strongly correlated metallic state is the Luttinger liq- 
uid (LL) in one-dimensional electron systems Q. The LL 
behavior has been observed in nanowires — of which the 
most prominent examples are the carbon nanotubes and 
the semiconductor quantum wires — through the power- 
law suppression 0, 3 of the tunneling conductance with 
decreasing bias voltage and/or temperature. 

When using the term "tunneling" , we often have at the 
back of our minds that the tunneling probes the proper- 
ties of the system into which the electron tunnels "nonin- 
vasively", i.e., the tunneling amplitude is infinitesimally 
small. A more subtle and complete understanding of the 
interplay of strong correlations and tunneling emerges 
when the latter is treated beyond the lowest order of per- 
turbation theory. It is the purpose of this paper to study 
the junction between a LL and a tunnel electrode (Fig.[T]) 
for arbitrary strength of tunneling. Apart from the con- 
ceptual interest, the "three-way junction" is considered 
a key element for device engineering in nanoelectronics. 

Our main result is the phase diagram for the 
interaction-induced renormalization of the parameters of 
the tunnel junction. We show that a homogeneous LL is 
actually generically unstable at low energies to arbitrar- 
ily weak tunneling and the breakup into two independent 
semi-infinite wires. The ZBA at the true (stable) fixed 
point (FP) is strongly enhanced. For sufficiently weak in- 
teraction or sufficiently strong tunnel coupling, the tun- 
nel conductance actually grows with decreasing energy 




FIG. 1: Junction between a quantum wire (horizontal) and a 
tunnel electrode (vertical). The arrows denote the scattering 
processes (and their amplitudes) for incoming electrons. 



scale (temperature, bias voltage) before it reaches maxi- 
mum and starts to renormalize toward zero. 

Let us specify the model. The Hamiltonian reads 
H = H,^ + He + Htnn, where (ft = 1) = 
Jdx (^—ifwipjj^dxipfj. + ^VoUfj^n-^) describes the LL 
wire (for compactness, we focus here on the spinless case), 
H — ± denotes electrons moving to the right (-I-) and to 
the left (— ). The electron operator -0, defined in the com- 
plete basis of scattering states of the three-way junction, 
is decomposed inside the wire as a sum of chiral compo- 
nents: ip{x) — V'/i(2^)j and n^[x) is the density fluc- 
tuation in the channel ^. We assume that the interaction 
potential between electrons in the wire is screened by a 
nearby metallic gate and take it to be point-like with the 
zero-momentum component Vq. The interaction between 
electrons with the same ji is then fully incorporated in the 
renormalization of the velocity v Q . We focus in this pa- 
per on the case of small a = Vq /2ttv. We assume that the 
interaction is present in a finite region of length L around 
the tunnel contact, which models a LL wire connected 
to noninteracting leads. For simplicity, we represent here 
the tunnel electrode as a noninteracting semi-infinite wire 

described by H^ = -iwc I],, M /°oo^2^ ^i^a^A" where 
are the chiral components of the electron operator ip{y) 
on the half-axis of y. In the absence of tunneling, the 
scattering states at y < are characterized by the phase 



FIG. 2: Renormalization group flows (schematic) for the point 
tunnel junction. Two stable FPs at p = and p — 1 describe 
the homogeneous wire and the breakup of the wire into two 
independent semi-infinite pieces, respectively. Tunneling into 
the wire is blocked at both points. 

of the reflection coefficient r — e*'^'' . 

The term i/tun describes the tunnel junction. We start 
by considering the simplest and commonly used model 
( "point contact" ) for the tunneling Hamiltonian: i?tun — 
t()4'Hy = = 0) + H.c, take to to be real and 

put the phase 4>r — 0. In the absence of interaction, 
the scattering amplitudes (Fig. [T]) can be shown to obey 
4 = -p, t,n = l- p, tout = t = -isgn(io)[2p(l - p)Y^^, 
r = 1 - 2p, where p = 2|to|VHc + 2|toP)- We see 
that the virtual transitions from the wire into the tunnel 
contact and back lead to backscattering in the wire. 

As is commonly known [Hj , a weak backscattering am- 
plitude in a LL is renormalized by interaction and be- 
haves as {A/\e\y~^ , where the energy e is counted from 
the Fermi level, A is the ultraviolet cutoff, and the Lut- 
tinger parameter K = {1 — ay^'^{l + a)~^l'^ ~ 1 — a 
for small a. At first glance, one might expect that the 
amplitude t\y is renormalized similarly. However, the fact 
that unitarity in the tunnel junction is imposed on the 
3x3 (and not 2x2) scattering matrix has dramatic con- 
sequences for the renormalization. We find (see below) 
that, to second order in a, the exact-in-p beta-function 
/3(p) = dpIdL, where £ = ln(A/|e|) for |e| > w/L, for 
the point tunnel contact reads: 

/3(p) =«p(l-p)[p-a(l-2p)(l-p + p2)/2] . (1) 

For p ^ 1 it reduces to 

/3(p) ~ -o?pI2 + ap\ p « 1 . (2) 

It is most important that /3(p) in Eq. ([2|) does not contain 
the term ap linear in both a and p, which would give 
the renormalization of similar to that for the weak 
impurity in a LL wire. To better understand this, recall 
that the renormalization can be described in terms 
of scattering off the Friedel oscillations produced by the 
scatterer. A subtlety of the point tunnel junction is that 
both tb and tin — in contrast to the impurity — are real 
and, as a result, the contributions to /3(p) at order 0{ap) 
from the Friedel oscillations that "dress" the junction at 
X < and x > exactly cancel each other. 

The first term in Eq. ^ means the usual suppression 
of the tunneling density of states oc (|e|/A)(i--^)'/2^ 
in a homogeneous LL. The second term comes from the 
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FIG. 3: Lowest-order scattering processes leading to the 
renormalization of the reflection amplitude r. Dots: tunnel- 
ing amplitude to. Wavy lines: bare interaction a between the 
right and left movers. The straight single and double lines de- 
note the electron Green function in the wire and in the tunnel 
electrode, respectively. 

renormalization of backscattering along the wire at first 
order in a, as found earlier in Ref. 7]. What seems to 
have not been discussed in the literature is that the two 
terms have opposite signs (for a > 0), so that (3{p) for 
small a vanishes at p = pc — a/2, which signifies an 
unstable FP (Fig. ^ leading to the phase transition sep- 
arating two phases with p = and p = 1 [l] . For p, a <C 1 
the solution of Eq. ([2]) gives 

Po 

2po/a+(l-2po/a)(A/le|)"V2 ' ^) 

where po is the bare value of p. 

We see that the tunneling is suppressed by interaction 
only if its bare amplitude is sufficiently small (or, equiva- 
lently, if the interaction is sufficiently strong). Otherwise, 
according to Eq. ([T]), the tunneling constant to mono- 
tonically increases in the course of renormalization. The 
tunneling transparency Gt — 2|t |^, however, shows a non- 
monotonic behavior with increasing to'. 

Gt = 4p(l - p) = Svv.ltoWiw, + 2\to\y . (4) 

If Pc < 1/2, Gt grows until it reaches maximum Gt — 
1 (which means a completely transparent contact with 
r = 0)atp=l/2 and decreases at larger p, eventually 
vanishing at p = 1. By contrast, Gt for p < Pc decreases 
monotonically in the course of renormalization. 

It is instructive to analyze the scattering processes dia- 
grammatically in the energy-space representation; in par- 
ticular, this makes it easy to count powers of to and a. 
The first and second terms in Eq. ([2]) correspond, respec- 
tively, to diagrams (b) and (a) in Fig.|3]for the amplitude 
r. In Fig. [3lja), tunneling with the amplitude to into the 
right-moving state is followed by scattering at first or- 
der in a off the Friedel oscillation, whose amplitude is 
proportional to tf, ~ 0{t^). Returning to the tunnel elec- 
trode costs one more power of to- Altogether, this gives 
ajtol^, which explains the origin of the second term in 
Eq. ([2]) (and — since the process of first order in a nec- 
essarily includes backscattering at the contact — it also 
explains why the term of order ap is absent). The other 
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FIG. 4: Skeleton diagrams (a) and (b) contribute to the exact- 
in-p /3-function at second order in a, diagrams (c)-(e) do not. 
The thick straight hnes denote the noninteracting electron 
Green function fully dressed by the tunneling vertices. 



scattering process [Fig. ^h)] is of second order in a be- 
cause of the creation of an electron-hole pair, but requires 
only two powers of to (only to get into the wire and come 
back): this gives a^jioP- 

More complicated diagrams of higher orders in p are 
constructed similarly: to second order in a they are com- 
pactly represented in Fig. |4j where the thick lines denote 
the noninteracting Green function dressed in all possi- 
ble ways by the tunnel vertices. Importantly, only one- 
loop diagrams (a) and (b) in Fig. |4] contribute to the 
/3-function to second order in a. Diagrams (c)-(e) cancel 
all singular interaction-induced terms in the ^-matrix ex- 
cept those given by the renormalization group equation. 
In particular, the perturbative corrections to r which 
come from Figs. 4(a) and (b) read 5r(^4a) — ^ —otttouttlC 
and(5r(4b) ^ -a'^Uout{tl\tb? +t*in\'ttn\'^)C/2. Calculating 
diagrams (a) and (b) for other amplitudes gives similar 
expressions, which — when written for the point tunnel 
contact in terms of p — all reduce to the single equation, 
Eq. (1). 

As follows from Eq. ([T]), tunneling into the wire is 
blocked at both FPs p — Q and p = 1: Gt vanishes as 
|e|" at p -> and as |e|" at p — > 1. We see that the 
ZBA for small a is strongly enhanced in the latter case. 
The difference in the exponents reflects the difference in 
the state of the wire: although the tunnel contact is de- 
coupled from the wire at both FPs, the wire at p = is 
homogeneous, while at p = 1 it is broken up into two dis- 
connected pieces (the exponent coincides then with that 
for tunneling into the end of a semi-infinite wire [2]). 

We now turn to a more general form of the tunnel cou- 
pling and show that the FP at p = 0, obtained above for 
the commonly used model of the point contact, is actu- 
ally generically unstable. The general parametrization of 
the 3x3 S'-matrix obeying time reversal symmetry and 
mirror symmetry x o —x (in particular, at an extended 
tunnel contact) gives for the moduli (which we focus on 
here) of the scattering amplitudes: 



\tb? - , W = l<ond' = 2p(cOS(^ - p) 



(5) 



where p and Lp are constrained by the condition < p < 
cos<y9 < 1 and the unitarity condition reads 2|tp + |rp = 
\tout\^ + l^inl^ + l^bP = 1- The point contact with real 
to corresponds to Lp = 0. The wire decoupled from the 
tunnel electrode corresponds to p = cos p. 



!5" 




FIG. 5: Renormalization group flows for the tunnel junction 
in the Gt-Gw plane for a — 0.15. The only stable FP is 
at Gt = Gw ~ (breakup of the junction into three parts). 
The unstable FP at Gt — Gt corresponds to the point pc in 
Fig. [2] The dashed line separates the regions in which Gt 
grows (above) and falls off (below) with decreasing energy. 



Diagrams (a) and (b) in Fig.|4l where the thick lines de- 
note now the "dressing" by the tunnel vertices described 
by Eqs. yield a set of two dosed flow equations 
for arbitrary p and tp. In the limit of small p ^ a, cos p, 
the /3-functions read ^d\: dp/dC ~ Q;p(sin^ p — ck/2) and 
dp/dC ~ {a/2) sin 2p. It follows that the FP (p = 0, (/3 = 
0) is actually a saddle point in (p, p) space. The key dif- 
ference from Eq. ([1]) is that the function dp/dC for small 
p now acquires the term apsin^t^. This means that at 
</? 7^ the cancellation of scattering on the Friedel oscil- 
lations at order 0{ap) [cf. the discussion below Eq. ^] 
is no longer exact. Written in terms of the dimensionless 
conductances Gt = 2\t\'^ and Gw — \tin\ 
flow equations are: 



*IV2 lUj, the 



dGt/dC = aGt{h+af2) 
dG^/dC = a(/3 + a/4) , 



(6) 



where /i = -l + Gt/2 + Gw, h = -l/2 + Gt(3-Gt)/8 + 
Gw(l ~ Gw), /a — Gt(l 4- Gw)/4 — 2Gw(l — Gw), fi — 
(1 - 2Gw)[Gt(4-hGt)/32-Gw(l - Gw)]. The flow stops 
at e given by T, bias, or v/ L, whichever is larger. 

The flow diagram for Gt and Gw is shown in Fig. [S] 
Two limiting curves are the line Gt = 4Gw(l — Gw) which 
describes the point tunnel junction, with the unstable FP 
from Fig. [2] at Gt = Gj: ~ 2a, and the line Gt = which 
describes a decoupled tunnel electrode. The FP (Gt = 
0, Gw = 1) is seen to be unstable to the decrease of Gw. 
It is this point that describes the commonly expected 
outcome of the renormalization: a homogeneous LL with 
tunneling blocked by the ZBA. In fact, this point has a 
finite basin of attraction (the line < Gt < Gj) only in 
the model of the point tunnel contact (which, albeit being 
widely used, is not generic in this sense). Generically, all 
flows are toward the point Gt = Gw = which describes 
the breakup of the junction into three disconnected parts. 

For 1 — Gw ^ 1 and the bare values Gto, 1 — Gwo ^ a. 
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FIG. 6: Left panel: schematics of the renormahzation of the 
tunnel conductance Gt on the log-log scale for the generic 
tunnel contact with 7^ (solid lines) and the point contact 
(dashed) . Curves I and IV (II and III) correspond to the bare 
conductance Gto below (above) the dashed line in Fig. [5] The 
arrow indicates the point at which the contact is reflection- 
less. The flow of Gt stops at the largest scale of e among T, 
v/L, and the voltage between the tunnel electrode and the 
wire. Right panel: the tunnel conductance of the point con- 
tact for a — 0.4 and Gto above (upper curve) and just below 
(lower) the critical conductance Gt — 0.45 (dashed line). The 
criticality retards the development of the ZBA. 

Eqs. ^ can be linearized, which gives: 

1 - Gw - Gt/4 = {Gto/GtT (1 - Gwo - Gto/4) , (7) 

where k ~ 2(2 — a) / a. We see that if the flow is initially 
close to the limiting curve 1 — Gw — Gt/4 predicted by 
the point contact model, it shows eventually a kink, af- 
ter which the "deviation" from the point contact model 
begins to grow sharply (k ^ for small a), see Fig. [S] 
The critical exponents change at the kink: before it both 
Gt and 1 — Gw decrease as |e|"^^^ [l2|> after it Gt con- 
tinues to follow this power law but 1 — Gw increases as 
|e|-2«. Near the FP at Gw = 0, Eqs. © predict that the 
exponent of Gt also changes and both Gt and Gw vanish 
at the FP as |e|" [H. The scaling of Gt,w oc |e|" de- 
scribes tunneling into the end of a semi-infinite LL wire 
[l3 | , which means enhancement (curve I in Fig. [H]) of the 
ZBA as compared to the conventional picture of tunnel- 
ing into a homogeneous LL (curve IV). If the bare values 
Gto, 1 ~ Gwo lie above the dashed line in Fig. [51 inter- 
actions first make the tunnel contact more transparent 
(curves II and III in Fig. (H]), so that both Gt and 1 — Gw 
first grow, but eventually the flow is attracted to the 
same FP Gt = Gw = 0. Note that the development of 
the ZBA can be strongly hindered in the vicinity of G^, 
as shown in Fig. [5] (right). 

To summarize, we have shown that the picture of tun- 
neling into a LL is qualitatively modified when the tun- 
neling amplitude is not treated as infinitesimally small. 
The conventional FP has a finite basin of attraction only 
in the model of the point tunnel contact, but taking a 
finite size of the contact (or any perturbation induced by 
the contact in the wire) into account makes it unstable. 



Generically, at the only stable FP the junction breaks 
up into three disconnected parts. Flowing toward this 
FP, the tunnel conductance may behave nonmonotoni- 
cally with bias voltage or temperature. Our predictions 
can be verified by systematically varying the distance to 
the tunnel electrode in experiments on carbon nanotubes 
or semiconductor nanowires. 
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